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ABSTRACT 


Computational  experiments  have  been  conducted  on  a 
digital  computer  to  illustrate  numerically  the  usefulness  of 
partial  coherence  techniques.  The  numerical  results  obtained 
demonstrate  the  necessity  of  employing  partial  coherence  functions 
when  a  multi-input-single  output  linear  system  is  involved. 

Two  seismic  noise  traces  were  randomly  selected  from 
available  Seismic  Data  Lab  sources  and  were  combined  in  various 
ways  to  obtain  two  correlated  input  traces  and  an  output  trace. 
From  this  data,  gain  factors  and  coherence  functions  were  com¬ 
puted  in  two  ways.  First,  just  one  input  and  one  output  were 
employed  in  the  computations.  Second,  both  inputs  were  ac¬ 
counted  for.  The  numerical  results  vividly  illustrate  the 
biased  answers  obtained  when  only  the  single  input  and  output 
are  considered.  This  is  contrasted  with  the  correct  results 
obtained  via  the  proper  model  of  a  two  input  system. 

These  procedures  were  also  extended  to  a  three  input- 
single  output  system  where  the  third  input  was  ignored.  In 
this  case  the  application  of  partial  coherence  functions  as 
an  indicator  of  the  existence  of  the  third  neglected  seismic 
data  source  is  ..1  lust  rated  . 


SEISMIC  PARTIAL  COHERENCY  STUDY 


1.  INTRODUCTION 

This  report  presents  the  results  of  a  preliminary  inves¬ 
tigation  into  the  application  of  partial  coherency  techr  ques 
to  the  problem  of  processing  seismic  data-  Ultimately,  these 
techniques  will  be  applied  to  the  detection  of  the  presence  or 
absence  of  unidentified  components  in  seismic  noise,  determina¬ 
tion  of  filter  responses  between  two  time  series,  and  description 
of  seismic  noise  fields. 

A  specific  application  of  partial  coherencies  to  seis¬ 
mology  would  be  the  modeling  of  a  seismic  noise  field.  If  a 
noise  field  has  only  one  component  propagating  across  an  area, 
then  the  output  of  one  seismometer  should  sufrice  to  predict  the 
output  of  a  second.  If  there  are  two  components  in  the  field, 
then  two  element  outputs  are  required  to  predict  a  third.  If 
there  are  n  components,  n  elements  are  required  to  predict  the 
output  of  one  additional  element.  As  will  be  seen  below,  the 
methods  of  partial  coherencies  can  be  used  to  determine  when 
sufficiently  many  inputs  are  being  used  to  predict  an  output. 
Thus,  it  will  be  possible,  by  using  partial  coherencies,  to 
determine  the  number  of  major  noise  components  present  and  the 
minimum  number  of  elements  required  to  study  the  noise  field. 

However,  before  (he  development  of  a  general  program,  it 
was  necessary  to  study  partial  coherence  techniques  applied  to 
the  solution  of  relatively  simple  problems  under  controlled 
condit ions  . 

The  application  of  power  spectra  and  cross-spectra  to 
determine  frequency  response  functions  for  simple  linear  systems, 
where  a  single  input  and  output  are  clearly  defined,  is  wall 
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established.  Application  of  these  ideas  to  more  complex  systems 
such  as  the  earth  where  there  are  many  input  and  output  points 
is  not  so  well  known.  Thus,  a  brief  review  of  these  topics  is 
presented  below. 

Consider  the  case  of  a  simple  linear  system  with  a  well- 
defined  single  input  (excitation)  point  and  single  output 
(response)  point.  If  the  input  x(t)  is  a  stationary  random 
process  with  zero  mean,  then  the  output  y(t)  is  also  a  station¬ 
ary  random  process  with  zero  mean.  Now,  if  G  (f)  and  G  (f)  are 

x  y 

the  one-sided  power  spectral  density  functions,  and  G^^(f)  the 
cross-power  spectral  density  for  the  input  and  output;  then  the 
frequency  response  function  for  the  linear  system,  H(f)  is 
completely  determined  by  the  two  relations 


G  (f)  =  i H ( f ) | 2  G  (f) 

y  x 


Gxy{f)  =  H(f)  Gx(f) 


(1) 

(2) 


A  quantity  of  particular  importance  in  more  complicate  i 
situations  is  the  coherence  function  which  is  defined  by 


G  (f) 
xy 


V  {  f  )  =  ■ 

'xy1  G  (f)  G  (f) 


(3) 


It  is  easily  shown  that  0  5  Y  (f)  5  1  ,  and  for  a  linear 

2  xy 

system,  y  (f)  =  1.  Thus,  the  coherence  function  may  be  thought 
of  as  a  measure  of  linear  relationship  in  the  sense  that  the 
function  attains  a  theoretical  maximum  of  unity  for  all  frequen¬ 
cies  in  a  linear  system.  Therefore,  if  the  coherence  function 
is  less  than  unity,  one  possible  cause  might  be  the  lack  of 
complete  linear  dependence  between  the  input  and  output. 
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Another  application  of  coherence  functions  to  single 
input-single  output  systems  is  to  determine  the  effect  of 
measurement  noise  on  frequency  response  function  measurements. 
Without  going  into  the  mathematical  details,  the  coherence 
function  in  this  case  is 


2 

V 

xy 


(f) 


1 


1.  +  3  .  +  3  +3,8 

i  o  1  o 


5  1 


(4) 


where  3^  i;  the  measurement  noise  to  signal  ratio  at  the  input 

and  8  is  the  corresponding  ratio  at  the  output.  Hence,  if 
o 

the  cohe  ence  function  is  found  to  be  significantly  less  than 
unity,  one  possible  cause  might  be  that  the  measurement  noise 
effects  are  not  negligible  and  must  be  taken  into  account  in 
determining  the  frequency  response  function. 

As  indicated  above,  coherence  functions  are  useful  aids 
in  the  analysis  of  single  input-single  output  linear  systems. 
However,  the  area  of  major  application  is  in  the  analysis  of 
multiple  input-multiple  output  linear  systems  which,  for 
example,  could  represent  the  response  of  the  earth  to  internal 
excitations  as  measured  at  several  different  points.  This 
would  be  the  case  when  an  array  of  seismic  instruments  is 
employed.  If  it  is  assumed  that  the  earth  has  a  linear  res¬ 
ponse,  and  that  measurement  noise  is  negligible,  then  a  low 
coherence  function  between  two  p-ints  will  serve  to  indicate 
the  presence  of  other  factors  which  contribute  to  the  output 

i 

but  are  not  being  considered. 

In  classical  statistics,  if  {x  }  is  a  sequence  of  random 

variables,  then  the  residual  random  variables  {Ax^j  are  defined 

as  the  result  of  subtracting  a  linear  least  squares  prediction 

of  x  from  x  .  Using  these  residuals,  partial  correlation 
i  l 

coefficients,  p.  can  be  computed.  The  notation  p. 


stands 


!i‘ 


for  the  correlation  between  a  and  x.  when  the  effect  of  x,  has 

1  3  k 

been  removed.  Partial  correlation  coefficients  are  useful  in 
determining  the  true  correlation  between  two  random  variables 
which  might  otherwise  be  obscured  by  the  effect  of  the  third 
random  variable. 

2 

By  analogy,  a  partial  coherence  function  y, .  ,  can  be 

i;pk 

defined  as  the  coherence  function  between  x^(t)  and  x_.(t)  after 

a  least  squares  prediction  of  x  (t)  has  been  subtracted  cut. 

K. 

It  should  be  understood  that  the  subscripts  i,  j,  and  k  may 
indicate  either  inputs  or  outputs.  The  utility  of  partial 
coherence  functions  is  similar  to  that  of  partial  correlation 
coefficients.  That  is,  a  high  ordinary  coherence  between  two 
variables  could  indicate  a  linear  relationship,  but  in  reality 
this  may  be  a  spurious  relation  due  to  correlation  of  one  of 
the  variables  with  a  third  variable.  If  the  partial  coherence 
function  is  calculated,  the  more  realistic  low  coherence  would 
be  uncovered.  Alternately,  the  opposite  effect  may  occur  where 
the  partial  coherence  will  be  larger  than  the  ordinary  coher¬ 
ence.  This  occurs,  for  example,  when  two  separate  inputs  both 
pass  through  linear  systems  and  make  up  the  single  output. 

Then  the  partial  coherence  between  either  one  of  the  inputs  and 
the  output  is  unity  due  to  the  linear  relations,  but  the  ordi¬ 
nary  coherence  will  be  less  than  unity.  This  follows  because 
each  input  contributes  to  the  output  and  this  fact  is  not 
accounted  for  in  the  computation  of  the  ordinary  coherence 
function  between  the  output  ard  a  single  output. 
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2,  SIMULATION  MODEL 


In  order  to  investigate  partial-  coherence  techniques 
under  con^iciled  conditions,  two  demonstration  cases  were  devel¬ 
oped,  These  are  described  briefly  below 

Case  1.  Given  x  (t),  x  (t),  and  y(t)  =  a  x  ( t )+a*x  ( t-T ) . 

1  Z  1  i  z  z 

Determine  a,  ,  a„  ,  and  T  in  the  case  where  x,  , 
12  1 

x  are  independent  and  in  the  case  where  x  ,  x 
are  correlated. 

Case  2.  Given  x^(t),  x2(t),  and 

y(t)  -  axx  j  ( t  )  +  a;?X2  ( t-T  )  +  a  x  ( t )  . 

Determine  a^  ,  a^  ,  T  ,  and  the  presence  of  x^ 
under  various  conditions  of  dependence  and 
independence  between  ,  x  ,  and  X3  * 

Detailed  procedures  for  generating  cases  1  and  2  are 
presented  in  Appendix  I  and  are  not  repeated  here.  Equations 
used  „n  calculating  expected  results  are  also  given  in  Appenaix  I. 
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3 .  ILLUSTRATIVE  EXAMPLE 


To  show  the  effect  of  partial  coherence  techniques  in  the 
analysis  of  a  linear  system,  one  example  will  not  be  discussed. 
This  example  plus  five  others  are  reviewed  in  detail  in 
Appendix  III . 

Consider  the  linear  system  sketched  below. 


z(t ) 


Mult iply 
By  0.4 


Multiply 
By  0.6 


Delay  byj _ | 


sec . 


Two  independent  seismic  sources,  x^(t)  and  z(t),  are  combined 
ro  form  a  new  process,  x  (t.),  such  that 

x2(t)  =  0.4x  (t)  +  0 . 6z ( t )  ( 

Secondly,  x2(t)  is  delayed  in  time  by  0.2  second  relative  to 
x^(t),  and  finally  combined  with  x^t)  to  form  y(t).  Thus 

y(t)  =  x1(t)  +  x  (t  -  0.2)  ( 


The  power  spec  ra  of  x^(t),  x  (t),  and  y(t)  are  shown  in 
Figure  1.  Because  of  the  way  x^t)  and  y(t)  were  formed,  the 
cross-spectra  are  given  by 


G,Jf)  "  0 . 4G  ,  (f) 


( 


(B) 


G  ff) 
lv 

G2y(f } 


1  +  0.4  cos  (0.4rrf)  -  j  sin  (0.4rrf)  jj  G^(f) 

(9) 

0'4Gll^f)  +  Lcos  (9*4^rf)  -  j  sin  (0.4"f)  G2  (f) 


If  the  transfer  functions  between  the  two  inputs  and  the  output, 
denoted  by  H^(f)  an<3  are  comPurfc^  (see  Appendix  I  for 

appropriate  equation)  without  taking  into  accoun :  the  correla¬ 
tion  between  x^(t)  and  x^(t)  completely  erroneous  results  are 
obtained,  as  shown  in  Figure  2.  It  should  be  noted  that  trans¬ 
fer  functions  are  complex  quantities  xn  general.  They  are 
denoted  here  by  gain  and  phase  variables 

H ( f )  =  i H ( f ) iej^(f)  (10) 

where  j H ( f ) I  is  the  gain  factor  and  #(f)  is  the  phase  factor. 

The  corresponding  true  transfer  functions  are  shown  in  Figure  3. 
These  are  computed  using  the  proper  multidimensional  formulas 
which  take  into  account  all  the  known  variables. 

The  ordinary  coherence  functions  are  plotted  in  Figure  4. 
If  used  by  themselves  to  infer  the  nature  of  the  system  being 
investigated,  one  would  tend  to  conclude  that  high  measurement 
noise  was  present  or  that  nonlinear  effects  were  present.  The 
true  nature  of  the  system  is  exhibited  when  the  partial  coher¬ 
encies  are  computed,  as  shown  in  Figure  5.  The  fact  that  all 
three  functions  are  equal  to  one  over  the  frequency  range  indi¬ 
cates  the  true  linear  relationships  which  exist  between 
x  ^  ( t )  ,  x2 (t ) ,  and  y(t ) . 

.xs  mentioned  before,  mere  detailed  evaluations  of  this 
and  five  other  cases  are  presented  in  Appendix  III. 


Measured  System  Response  Functions 


True  System  Response  Functions 


4.  CONCLUSIONS  AND  EXTENSIONS 


The  demonstration  cases  investigated  through  the  aid  of 
computer  simulation  methods  illustrate  the  usefulness  of  partia1 
coherency  techniques  in  analyzing  multiple  input-multiple  output 
systems.  The  fact  that  correct  values  of  transfer  functions  and 
coherence  functions  are  obtained  when  all  inputs  are  properly 
taken  into  account  has  been  numerically  illustrated. 

The  next  step  in  the  investigation  will  be  to  use  the 
computational  procedures  described  in  Appendix  II  to  compute 
transfer  and  coherence  functions  completely  from  actual  data. 
These  can  then  be  compared  with  the  expected  results  discussed 
in  this  report.  Questions  of  statistical  accuracy  which  will 
be  limited  by  sampling  errors  in  the  spectral  computations  can 
then  be  explored.  Another  area  for  future  stud’  is  the  trade¬ 
offs  between  statistical  accuracy,  number  of  inputs,  and 
computation  t’me. 

Finally,  a  general  program  must  be  developed  for  the 
multi -input  case  and  applied  to  the  measurement  and  interpreta¬ 
tion  of  seismic  noise  in  large  arrays. 
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EXPECTED  RESULTS  FOR  SEISMIC  COHERENCY  EXPERIMENT 


1 .  INTRODUCTION 

The  following  equations  and  input  data  form  the  basis  from 
which  expected  experimental  results  were  obtained  by  machine 
computation.  The  results  of  these  computations  provided  mathe¬ 
matical  checks  on  the  computer  simulation  work.  The  defining 
equations  for  Case  I  are  presented  in  Section  2  while  the  modifi¬ 
cations  required  to  run  Case  II  are  given  in  Section  3. 

The  two  demonstration  cases  were  based  on  data  constructed 
in  the  following  manner: 

I.  A  signal  of  seismic  noise  y(t)  is  composed  of  two  other 

known  statistically  independent  seismic  traces  x^(t)  and 

x„(t)  where 
£ 

y ( t )  =  a^x^ ( t )  +  a2X2't  “  ^ 

(a)  The  two  signals  x  (t)  and  x^ft)  are  independent. 

(b)  The  two  signals  x  (t)  and  x2(t)  are  statistically 
correlated. 

II.  The  signal  of  seismic  noise  y(t)  generated  as  indicated 
above  in  (I)  was  modified  by  the  addition  of  a  third 
component  x3(t).  For  this  case,  the  formula  for  y(t)  is 

y  ( t )  *  a  x1(t)  +  a2X2^t  “  +  a^U)  (2) 

(a)  The  three  signals,  x^t),  x2(t),  and  x^(t)  are 
independent . 

(b)  The  two  signals  x  (t)  and  x2(t)  are  correlatec*'  but 
x  3  ( t )  is  independent  of  x^t)  and  x2(t). 


I 
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(c)  The  two  signals  x^(t)  and  x^(t)  are  ^dependent, 
but  x^(t)  is  correlated  with 

(d)  The  two  signals  x^(t)  and  x^Ct)  are  correlated  and 
x^(t)  is  correlated  with  x^(t)  but  independent  of 


2.  CALCULATION  FOR  CASE  1 

A  block  diagram  of  Case  I  is  shown  in  Figure  1.  Referring 
to  the  diagram,  x^(t)  and  z^(t)  are  samples  from  independent 
random  processes.  The  process  x?( t)  is  formed  by  a  linear  com¬ 
bination  of  and  z^Ct),  thus 

x^Ct)  =  ax^Ct)  +  (1  -  a)z1(t)  (3) 

The  random  process  y(t)  is  obtained  by  applying  gain  factors  a^ 
and  a^  to  x^(t)  and  x^(t),  ''espect ively ,  and  by  delaying  x2{t) 
by  an  amount  T  so  that 

y  (t )  =  a^Ct)  +  a2x9(t  -  T)  (4) 


3.  • 

x}(t) 

and  z  ( t  >  independent 

b. 

x  2  ( t ) 

■  ax  (t)  +  (1  -  a)z  (t) 

c  . 

y(t)  ; 

=  a-jX  (t )  +  a2x2  (t  -  T) 

d  . 

H1(f) 

=  al 

*9  - 

H  2  (  f ) 

-]2nfT 

=  a  e 

2 

Figure  1. 

Slock 

Diagram  of  Method  for  Construction  of  Case  1  Signals 

For  computation  of  the  various  coherence  and  frequency 

response  functions,  the  following  format  was  used.  Notation 

used  follows  that  employed  in  Ref.  [lj. 

In  che  computations  described  below,  "true"  and  “hieasured" 

frequency  response  functions  are  calculated.  The  measured 

frequency  response  functions  H,  (f),  i  =  1,  2  are  defined  as 

lM 

the  functions  which  would  be  calculated  only  from  knowledge  of 
one  input  and  one  output.  These  functions  are  incorrect  when 
additional  inputs  exist  which  are  not  accounted  for .  The  true 
function  1L  (f),  i  =  1 ,  2  are  defined  as  the  versions  computed 
using  the  proper  multidimensional  formulas  which  account  for  all 
the  known  variables.  Hence,  the  difference  between  the  true  and 
measured  frequency  response  functions  will  reflect  the  advantage 
in  using  the  proper  computational  formulas. 


1.  Specify  G..(f)  and  G  (f)  in  the  frequency  interval 
11  zz 

0.5  to  2.0  cps  at  0.1  cps  increments. 


2  . 


4. 


5. 


Compute 


Compute 
G  (f) 

yy 

Compute 


G22(f)  using 

G__(f)  *  a2G  (f )  +  (1  -  a) 2G  (f)  (5) 

22  11  zz 

G  (f)  using 

yy  ( 6 ) 

[a^  +  2axa2a  cos  ( 2trfT )  J  G^f)  +  a^G^U) 

G12(f)  using 

G12(f)  *  aG11(f)  (7) 


Compute  G,  (f) 

iy 

C.  (  f ) 

iy 

Q,  (f) 

I  \7 


-  j  *  ,  (  f  ) 

*  C.  (f)  -  jQ.  (f)  =  Ig.  (f ) je  iy 
ly  ly  Iv 

=  La1  +  aa2  cos  ( 2nfT )  JG^  ( f )  (8a) 

=  aa2  sir  (?mfT)  g  (f)  (8b) 


I 


g  (f 5 i  - 

iy 


/  2  2 

'c .  (f)  +  Q,  (f) 

ly  iy 


(9a) 


i  Q,  (f) 

V£)  ■  tan  EfuT 

ly 


(9b) 


6.  Compute  G„  (f)  using 

2y 


(10a) 


C2y(f)  -  a1aG1]_  ( f )  +  a2  cos  (Jr£T)G22(f) 


Q2y ( f )  =  a2  sin  (2nfT)G22(f ) 


(10b) 


lG2y(f)l  %/C2y  <£>  +  °2y  (£i 


(11a) 


i  Q 9 

*2y(£)  *  tan'  ^fTf) 

2y 


(lib) 


7.  Compute  the  "true"  frequency  response  functions, 

HlT ( f )  and  H  T(f),  using 

G  ( f )C  (f)  -  G  (f)C  (f) 

ReH.f)^ - - 12  -  2y  U2a) 

Gn(£)G2,(f)  -  G12  (f) 


G,,(f)0,  (£)  -  G  ( f ) Q  (£) 

ImH  (  f )  -  - & - 22  2  lY  (12b) 

GU(£)G22!f)  '  G12  (£) 


G,,(f)C,  (£)  -  G.if)C(f) 

ReH  (f)  =  — ^ ^ - 12  .  ly--  (13a) 

Gu(f)G22(f)  -  G^  (£) 


G,  ( f )Q1  (£)  -  G  (£)Q  (£) 

ImH  (f)  -  -12 - lx- - ;1  2-2-y—  U3b) 

GU(£)G22(£)  -  G12  (£) 


1-5 


2y  V  ;  G  ( f ) G  (f) 

22  yy 

12.  Compute  the  following  residual  spectra 
GU.2(£)  “  -  ''122,£)] 

Gll-y(f)  *  Gll(£)[1  -  Y122<£>] 

G22-l'£)  =  °22(£,[1  -  V2y2(£>] 

G22-y(f)  =  G22<£>[1  -  v2y2  (  £)  ] 

G  .  (f)  =  G  ( f ) f 1  -  Y„  2(f) 

yy  i  yy  L  2y  J 

g  ,(f)  =  g  (f)ri  -  y  2(f)  i 
yy*2  yy  l  2y  J 


(20) 

(21) 

(22) 

(23) 

(24) 

(25) 

(26) 


13.  Compute  the  following  quantities:  (Note  that  the 

argument,  f,  has  been  dropped  from  all  functions  for 
simplicity . ) 


I 


7 


3.  COMPUTATIONAL  PROCEDURE  MODIFICATIONS  FOR  CASE  II 


The  block  diagram  for  Case  II  is  shown  in  Figure  2. 
addition  of  a  third  process,  z^lt)  which  is  combined  with 
to  form  x^(t)  is  the  only  change  from  Case  I. 

The  computational  procedure  presented  previously  was 
fied  as  follows  for  Case  II: 

2.  Compute  G22(f)  and  using 

G  (f)  =  a2G  n (f )  +  (1  -  a) 2G  (f ) 

22  H  Z1Z1 


G,,(£!  =  B2G,, (f)  +  (1  -  B)2G  (£) 
33  11  z2z2 


3.  Compute  G  (f)  usinq 

yy 


G  (f)  =  a  2  +  2a  a.  a  cos  (2nfT)  +  2a  a  6 

yy  Li  i2  i  j 


+  2a2a^  aS  cor  ( 2nfT )  JG^  (  f ) 


+a22G22(f>  +  a32G33(f> 


5.  Compute  G  (f)  =  C  (f)  -  ]Q  (f) 
ly  ly  ly 


C  (f)  ~  a  +  a  a  cos  (2rfT)  +  a  3  jG  .(f) 
ly  L  1  2  3  J  11 


Q  (f)  =  a  a  sin  (2nfT)G  (f) 
ly  2  il 


T  -  9 


The 
x1  (t ) 

modi- 

( 5 '  a ) 

(  5  *  b ) 


(6'  } 

(8  *  a) 

( 8  '  b ) 


y  (t ) 


a . 

x1(t),  z1(t),  z2(t) 

INDEPENDENT 

b. 

x2(t)  -  ax1(t)  +  (1  - 

a)  zx(t) 

c . 

x3(t)  =  3x1(t)  +  (1  - 

3)  z2(t) 

d. 

H  (f)  -  a 

e . 

h2(£,  -  ,2 

f . 

H3(f)  =  a3 

Figure  2.  Block  Diagram  for  Method 

of  Construction  of  Case  II  Signals 

w 


and 


iy 


iy 


c.  (f) 
iy 


6.  Compute  G  (f)  using 


(9‘a) 


(9'b) 


C2y(f)  =  +  a3aej  Gu(f)  +  a2cos(2nfT)  G^f)  (i^'a) 


Q2^(f)  =  a2sin  (2nfT)  G22(f) 


and 


(10'b) 


( 11 ' a ) 


(ll'b) 


All  other  computations  were  unchanged 


REFERENCES  -  APPENDIX  I 

1.  Enochson,  L.  D.,  "Frequency  Response  Functions  and  Coherence 
Functions  for  Multiple  Input  Linear  Systems,"  NASA  Cr-32, 
prepared  under  contract  NAS-5-4590.  April  >1. 


APPENDIX  II 

COMPUTATIONAL  PROCEDURES 
FOR  SEISMIC  COHERENCY  EXPERIMENTS 


m 


computational  procedures 

FOR  SEISMIC  COHERENCY  EXPERIMENTS 

The  computational  procedures  here  are  intended  to  give 
correct  results  when  the  proper  physical  model  is  a  two  input- 
single  output  time  invariant  linear  system.  This  model  is 
depicted  in  Figure  1. 


Figure  1.  Two  Input  Single  Output  Linear  System 

The  two  input  processes  x^(t)  and  x2(t)  both  are  assumed  to  be 
stationary  with  zero  mean  but  need  not  be  statistically  inde¬ 
pendent  . 

The  procedures  described  naturally  separate  into  three 
basic  steps.  There  are: 

(1)  Compute  auto  and  cross-correlation  functions. 

(2)  Transform  the  weighted  correlation  function  to 
obtain  the  spectral  densities  where  either  or 
both  the  Hanning  and  Parzen  weighting  functions 
have  been  applied. 

(3)  Compute  the  frequency  re  ,ponse  functions, 
ordinary  coherence  functions,  and  partial 
coherence  functions. 


II 


1 


The  discrete  digi  ized  points  of  the  zero  mean 


stationary  processes 

are  dt  ,oted 

as  follows: 

x1  (t ) 

: 

(i) 

i  =  1,2, - N 

X2!t) 

: 

x2  (  i) 

i  =  1 . 2  ,  .  .  .  ,N 

y(t ) 

* 

y(i) 

i  =  1 , 2  ,  .  .  . ,  N 

The  quantity  N  represents  the  number  of  discrete 

data  points  for  each  seismic  data  record.  The  symbol  m  will 

be  used  to  designate  the  maximum  number  of  lags  at  which  the 

correlation  functions  are  computed  lefined  by  m^T  =  t  .It 

max 

is  expected  that  the  digitizing  rate  is  20  cps ,  and  the  desired 

bandwidth  B  of  the  spectral  density  estimates  is  0.1  cps. 

Hence,  since  B  =  1/t  ,  one  must  have  m  ■*  200.  Also, 

max 

t  =10  sec  so  to  satisfy  a  requirement  of  10t  <  T,  a 

max  max— 

record  length  of  T  =  100  sec  would  be  required. 

1.  COMPUTATION  OF  THE  CORRELATION  FUNCTIONS  AND 

RAW  SPECTRAL  DENSITIES 

1 . 1  Autocorrelation  and  Cross-Correlation  functions 

N-t 

R11(T)  ’  N  l  X1(i)X1(i+T) 

i=l 

N-t 

R22<t)  x2(i)X2(i+T) 

i=l 

N-t 

R  (T)  ■  ”■  Y  y(i)y(i+T) 
yy  N 

i  =  l 


II 


2 


R12(T) 


N- ' 

I  7 

N  L 


x(i)x( i+T ) 


i=l 


R 


21 


(t) 


N-T 

1  v  ,  .  v 

-  .  x  ( i ) x  (i+T) 

N  /■  2  1 

i=l 


N  r 

R,  (t)  e  ~  )  x  (i)y(i+T) 
j.y  N  t-i  l 

1=1 


N-T 

R  . (t)  =  ^  )  y(i)x  (i+T) 
yl  N  '  i  1 

i=l 


R 


2y 


(T) 


N-T 

nZ  x2(i)y(i+T) 
i=l 


N-T 


R 


v2 


(  T) 


1  V 

N  Zj 


y(i)x  (i+T) 


i  =  l 

i  =  0,1,2,... ,m=200 

1 , 2  Determination  of  a  Peak  in  the  Cross-Correlation  Function 

For  each  of  R12(t),  R21^'  Rly^  '  Ryl^T^'  R2y^T^' 
and  R  „(t),  the  maximum  statistically  significant  peak  (if  one 

y2 

exists)  is  to  be  found  and  the  cron^-correlat ion  functions 
translated  appropriately.  This  search  for  a  significant  peak 
should  be  limited  to  an  interval  of  time^-T,Tj  which  corresponds 
to  maximum  physically  meaningful  time  delays.  This  translation 


II 


3 


of  the  correlation  function  corresponds  to  a  crude  "prewhiteninq* 
of  the  cross-spectral  density.  Hence,  if  some  more  general 
prewhitening  for  the  cross-spectral  computations  was  employed, 
this  translation  need  not  be  performed. 

Let  R^_.(t)  and  R_,^(t)  represent  the  two  parts  of  any  of 
the  three  cross-correlation  functions.  Note  that  R  ( -t )  * 

i  j 

R u  T  )  *  Determine  the  maximum  function  value  of  these  two 
parts.  That  is, 

"max  *  max  lRji<T  > I  '  T  *  O'1 . »]  <«> 

T  J  J 


Next  calculate  the  maximum  correlation  coefficient 

R 

r  =  _ _ 

R . . (0) 

D  3 


As  a  reasonable  test  of  statistical  significance,  it  is  suggested 
that  r  be  compared  against  P  =  0.25.  To  perform  this  test, 

ITloX 

compute 


max 


which  may  be  assumed  to  be  norma  ly  distributed  with  mean  £ 
and  variance  l/(n-3)  where 


II 


4 


1 

2 


1  +  0.25 


.255 


j 


n 


1  -  0.2 


and  the  effective  sample  size  is  given  by 

n  =  2BT 


It  is  assumed  that  the  record  length  T  will  be  about  100  sec 
or  longer.  The  signal  bandwidth  is  B  and  a  reasonable  value 
for  B  should  be  1.5  cps  in  which  case  r.  =  3T .  Hence,  if 


Z  =  tan  h”1  | r  j  >  .255  +  2.33 

max  —  — 

/3(T  -  1} 


.255  +  1 -34  (8) 


then  the  hypothesis  that  r  is  larger  than  0.25  or  a  negative 

max 

r  is  less  than  -0.25  is  accepted  at  the  1%  level  of  signifi- 
max  r 

cance .  The  quantity  2.33  is  the  99th  percentile  of  the  normal 
distribution  function. 


1 . 3  Translation  of  the  Cross-Correlation  Function 

Let  Tq  *  sq^t  represent  the  time  point  at  which  a 

statistically  significant  r  for  x,(t)  and  x  (t)  occurs. 

max  1  2 

This  value  is  t  =  0  if  no  significant  peak  in  the  sample 
cross-correlation  function  exists.  Then,  the  quantities 
A12^T^'  B12^T^  35  ^e^ine<^  helow  are  computed: 


A12(T)  ‘  KR12(T  +  S0>  +  R21(T 


-  s 


1 


B  (  T  ) 
12  V 


KR12(T 


s  )  - 
0 


Rn( 


V] 


(9a) 

(9b) 


II 


5 


and 


represent  delay  values  determined 


i-et  '  =  s,ut 

L  X 


t  s  o  t 

2  2 


for  the  cross -correlations  between  x,(t)  and  y(t)  and  x_(t) 

1  2 

and  y(t),  respectively.  Also  compute 

Aly(T)  “  I  L*ly,T  +  V  4 

Bly(T|  -  I  lRly(T  +  sl>  - 


2y 


2y 


1  r 

—  R„  ( T  + 

2  l  2y  2 

T"  R„  (t  +  s_ 
2  l  2y  2 


"s 

yl'T  -  VJ 

( 1  Oa ) 

(10b) 

(  T  -  S  )  1 

y2  2  •  j 

(11a) 

y2lT  -  »2>] 

i  lib) 

2.  COMPUTATION  OF  WEIGHTED  FOURIER  TRANSFORMS 

Two  different  weighting  functions  are  to  be  allowed 
in  the  computation  of  the  power  spectral  density  function. 

These  are  the  Hanning  weighting  D^Ct)  and  the  Farzen  weighting 
D  (t).  This  weighting  may  be  performed  in  a  slightly  more  compu¬ 
tationally  efficient  method  by  smoothing  of  the  raw  spectrum 
estimates.  However,  for  convenience  in  varying  weighting 
functions,  it  would  seem  desirable  to  multiply  the  correlation 
functions  by  the  weighting  function. 

Hence,  compute  the  weighting  function  values 


D1(t  ) 


1 

2 


cos 


t=0, 1 , 2 , . . .m  ( 12 ) 


II 


6 


r  -i 


v. 


t  =  0,1,2,... ,m/2 


(13) 
,  m 


Compute  the  spectral  density  estimates  by  Fourier 
transforming  the  correlation  functions.  If  the  necessary 
sine  and  cosine  values  for  the  transform  are  not  precomputed 
and  stored,  they  should  be  generated  from  the  following 
recursion  formula  where  sin  (n/m)  and  cos  (n/m)  are  necessary 
to  get  started. 


”cos  (i  + 

1 )  n/rrf 

cos  in/m  -sin  in/m 

cos  n/m 

^sin  (i  + 

1 )  tt/itl 

_sin  in/m  cos  in/m. 

„s  i  n  n/m_ 

2 

l—l .  2  _ 

The  Fourier  transform  equations  for  the  autocorrelation  functions 
are 

m-1 

r  TTTk  1 

G  (k)  =  2ut;r(0)  +  2  )  D  (t)R(t)  cos  -  .  k=4 , 6 , - 34 

1  u  L.  m-i 

T  =  1 

where  it  is  tacitly  understood  that  R(t)  stands  for  any  of  the 
three  autocorrelation  functions  R^(t),  R^2^t^'  and  Ryy^T^' 

The  evaluation  of  the  above  equation  is  therefore  performed 
a  ’'~tal  of  six  times?  two  different  weighting  functions  for 
each  of  three  autocorrelation  functions.  The  ^oectral  density 
is  evaluated  at  the  16  values  of  k  which  correspond  to  frequencies 
f  =  0/2,  0.3, - 1.7. 

The  transform  for  the  cross-correlation  function  consists 
of  a  real  part  c(k)  and  an  imaginary  part  Q(k).  The  equations 


i 


ii 


7 


employ  the  special  A  s  and  B's  defined  by  Eqs  ,  (9),  (10) , 

and  (11).  If  G  (k)  represents  the  cross-spectral  density 
between  some  arbitrary  u  and  v,  then 


Re  ;  G  { k)  1  =  C  (k/  =  2di 
L  uv  j  i 


m-  1 

c’  tttk 

A (  0 )  2  AiTjD  (t)cos  - - 

i~>  i  m 

T=1 


(16) 


Im  G  ( k  J  -  Q  ( k 
L  UV  J  1 


m-1 

V  TTTk 

4di  )  B  i  t  j  D  (  T  j  s  m  - - 

L  x  m 

i  A 


k  =  4,  6 , 


..  34 


Again,  Eqs.,  (16)  are  evaluated  a  total  of  six  times.  The 

quantities  A(i)  and  B(t>  represent  an<^  B  ^  2  iT  1  '  * 7  ) 

and  B,  (t)  and  A„  (tj  and  (t).  These  quantities  are  each 
ly  2y  2y 


used  twice,  once  with  D. 


and  once  with  D (  T  , 


1  2 
At  this  point .  all  the  necessary  spectral  densities 

have  been  evaluated.  These  are 


G 


11 


k) 


G  Ck 

yy 


G12U) 


12 


;a!2'k 


G.  (k)  =  C  (k)  -  jQ,  \k) 
ly  ly  ly 

G  Ck)  =  C_  (k)  -  jQ  (k 
2y  2y  2y 

Each  of  the  quantities  e^e  available  in  two  forms  corresponding 
to  the  two  different  weighting  functions. 

At  this  point.,  the  cross-spectral  densities  should  be 
adjusted  to  account  for  the  time  translation  that  was  applied 


II 


8 


to  the  cross-correlation  functions.  This  translation  by 


t  ime 
by 


amounts  to  a  multiplication  of  the  cross-spectrum 


j2nfs^T 

e 


Hence,  the  sample  cross-spectrum  must  be  multiplied  by 


-  j2rrfs  .  ut 

i 

e  to  readjust  .  The  adjusted  cross-spectra  will  be 


temporarily  denoted  with 

a  prime  ( ‘ ) * 

The  equations  are 

1 

C 

'l2(k) 

=C12(k) 

cos 

2  tt  f  s  a  T  - 
o 

Q12(k 

si  n 

2-£s0At 

1 

Q 

'l2(k) 

=C!2(k) 

s  in 

2nfs^^T  + 

Q12tki 

cos 

C 

ly(k) 

=  ciy(k) 

cos 

2-fs^T  - 

Q.  (k) 
iy 

sin 

2TTf31tiT 

Q 

iy<k) 

=  c  (k) 
iy 

sin 

2^3^^  + 

Q,  (k) 
iy 

cos 

2rrf  s^At 

! 

C 

( k) 
2y 

=  C_  (k) 

2y 

cos 

2~fs  At  - 

Q_  (k) 
2' 

sm 

2~f 

Q 

'  (k) 

2y 

=  C.  (k) 

2y 

s  in 

2TTfs.AT  + 

Q-.  (k) 
2y 

cos 

2TTfS2£>-i 

Note  that 

these 

computat ions 

are  performed  at 

the 

usual  16 

frequency  values.  Subsequently  the  primes  will  be  dropped 

and  the  adjusted  spectra  will  simply  be  denoted  by  G  (k)  r 

G.  (k) ;  and  G  (k) . 
ly  2y 


II 
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COMPUTATION  OF  FREQUENCY  RESPONSE  FUNCTIONS ,  ORDINARY 
COHERENCE  FUICTIQN^  ,  AND  PARTIAL  COHERENCE  FUNCTIONS 


The  computations  to  follow  will  y l"'  a  "true”  and 

"measured”  frequency  response  function  H  (k)  and  H  (k)  as 

IT  iM 

defined  in  the  discussion  on  experimental  results.  Then  the 
ordinary  coherence  function  formulas  will  be  given  followed 
by  the  partial  coherence  function  formulas.  It  is  to  be 
understood  that  all  formulae  are  evaluated  twice,  once  for 
each  weighting  function. 

Several  of  the  subsequent  formulas  correspond  exactly 
to  formulas  appearing  in  the  write-up  on  expected  results 
calculations.  In  each  case  where  this  occurs  a  special 
mention  is  made  so  that  separate  subroutines  could  be  common 
to  both  programs . 

3 . 1  Frequency  Response  Functions 

(1)  First  compute  the  "measured"  frequency  response- 
functions.  These  formulas  correspond  to  Step  8  of  the  expected 
results  computations. 

H(k)  =  ReH(k)  t  j  In  H(k)  =  |h(3c)  | 


k  =  4,  6 . 34  (18) 


ReVk) 


Clv(k> 

Gu(k) 


Q.  (k) 

iv 

Gll(k) 


II  -  10 


. . . """ . . . . . '  . ! . . . 


Compute  the  gain  and  pha.se  factors 


|H.M(k)|=^/rReH  (k)]2  +  ^ImH  (k) 

1  v  L  1M  J  i_  1M  J  iM 


-1 

ReHlM!k) 


In  subsequent  formulas  the  k  will  be  dropped  for  notational 
convenience  and  it  will  be  understood  that  all  computations 
are  performed  at  16  values,  k  =  4. 6,... ,34. 


{ 19) 


Compute  the  gain  and  phase  factors 


/ - ; - - - -  1  ImH2M 

HlM(k)|  -V(ReH2M)  *  (Mal  .  'tan  j— - 


(2)  Computation  of  the  "true’'  frequency  response 
function.  First,  compute  the  squared  magnitudes  of  the 
cross-spectra  for  later  use. 

■  .2  2  2 


G12' 

=  C12 

+  Q 

12 

Giy'2 

2 

ciy 

2 

+  Qiy 

(20) 

%'2 

2 

C2y 

2y 

II 
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(4)  Computation  of  the  partial  coherence  functions 
First,  the  residual  spectra  are  computed  which  corresponds 
to  Step  12  of  the  expected  results  computations. 


G,  ,  •  _  “G,.,  1  -  v ,  . 


11  •-  2  11  L 


12  J 


11  ‘  y  11 


G,  ,  r  1  -  v.  2 


°22' 1  °22  [  1  V12  J 


1  -  v. 


G  *  G 

22‘ y  22 


[  1  -  v2y  J 


(25) 


G  -  G  1  -  V. 

yy  i  yy  l  iy  j 


G  —  G  1  -  V 

yy  2  yy  L 


2y  j 


Compute  the  residual  cross-spectra: 


G-.  t  =  t  "  jQ-,  T 

12  y  12  y  12  y 


C 12 ‘ y  C 12 


C.  +  Q. 

ly  2 y  ly  2’ 


(26) 


Q12’y  "  °12 
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Compute  the  partial  coherence  functions.  This  is  identical 
to  Steps  14,  15,  and  16  of  the  expected  results  computations 


Y12'y 


G12  ‘  v 
Gll*  y°22 ' y 
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^*2  “  Gll-2Gyy2 


(29) 
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This  concludes  the  basic  computations.  The  output  format 
duplicates  that  of  the  expected  results  computations.  This 
allows  a  one-to-one  comparison  of  the  experimental  results 
with  the  theoretically  expected  results.  In  addition,  the 
power  spectral  density  functions  G,,(f),  and  G  (f) 

ii  2  22  2  yy 

and  the  partial  coherence  functions  (f),  v,  „(f), 

2  12' y  ly  2 

V*-  .(f)  are  output  to  comp  re  with  known  values. 

2y '  1 


evaluate  nes  and 
cosines  fox  trans¬ 
form  recursively 
using  Eq.  (14) 


evaluate  weighting 
function  D.(r)-see 
Eqs.  (12)  and  (13) 


compute  weighted 

cosine  transforms  of 

three  autocorrelation 

functions  to  obtain 

(k) ,  G  (k) ,  G  (k)  , 
22  yy 

using  Eq .  ( 15 )  . 


compute  cosine  and  sine 
transforms  of  cross¬ 
correlation  functions  to 
obtain  real  and  imaginary 
parts  oi  cross-spectra 
ci2(k),  cly(k).  c2y(k), 

G10(k),  Q  (k) ,  and  Q  (k) , 
12  ly  2y 

using  Eqs .  ( 16 ) . 
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Readjust  cross-spectra  for 
time  translation.  Compute 


C12(k) 

Q12(K) 

Cly(k) 

Q|y(k) 

c’  (k) 

2y 

Q0  (k) 
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C12(k) 

Q12(k) 
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Apply  Eqs  .  { 1"?  >  . 


Compute  ReH  .  ImH  ,  , 

lM  }M 

^  H  jM  ^  '  *  jM 

"measured"  frequency 
response,  j  =  1,  2. 
Use  Eqs.  (18) ,  (19) . 


Compute  |G 


12 


|G 


lyl  '  lG2y^  '  and  *' 


Use  Eqs .  (20) ,  (21) 


Compute  ReH.m,  ImH,  , 
jT  jT 


H  ,  $  .  -  "true" 

jT  jT 

frequency  response,  j  =1,2 
Use  Eqs.  (22) ,  (23) . 
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compute  ordinary  coherence 

2  2  2 

functions  v, _ ,  v,  ,  and  v„ 
12  1 y  2  y 

f  rom  Eqs  .  ( 24 )  . 


compute  residual  power  spectra 

G  ,  G  ,  G  .  G  , 

11*2  11*  y  22*1  22  3  y 

G  , . ,  and  G  . _  from  Eqs .  (25). 

yy  1  yy  2 


compute  residual  co  and  quad- 
spectra  C12>yf  Q12.y.  C^, 

Qly*2'  and  C2y’ 1 '  Q2y* 1  fr°™ 
Eqs .  (26)  ,  (27  )  ,  and  (28)  . 


compute  partial  coherence 

^  2  2 
function  Y,  ,,  »  V.  ,  , 

12  y  ly  2 

and  y~  . using  Eqs.  (28). 
2v  1 


prepare  for 

printer  output  tape: 

G,  ,  , 
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11 

22  yy 
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jM 
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dt 
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printer) 
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Is  this  case  to  be  plotted 
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APPENDIX  III 

ANALYSIS  OF  RESULTS  OF  COMPUTATIONAL  EXPERIMENTS 


analysis  of  results  of  computational  experiments 


1 .  INTRODUCTION 


The  purpose  and  general  results  of  the  computational 
experiments  are  described  in  the  main  body  of  the  text,  but 
certain  pertinent  details  will  be  reviewed  here.  The  over-all 
objective  of  the  computational  experiments  is  to  develop  a 
numerical  feeling  for  the  magnitude  effects  of  various  data 
parameters  on  the  results  of  interest.  A  side  benefit  is  gaming 
familiarity  and  experience  with  the  general  computational 
procedures  required  in  the  analysis  of  multi-dimensional  linear 
systems . 

These  preliminary  numerical  experiments  provide  results 
which  are  not  influenced  by  statistical  variability  so  that  an 
over-all  impression  of  the  effectiveness  of  the  theory  may  be 
gained.  The  computations  are  based  on  actual  seismic  data  so 
that  some  of  the  physical  details  of  the  problem  are  reflected 
in  the  experiments  even  though  the  statistical  variability 
which  is  encountered  in  actual  data  analysis  has  been  avoided 
in  these  initial  evaluations. 

As  mentioned  in  the  main  body  of  the  report,  two  basic 
cases  were  examined.  Case  I,  a  two  input-single  output  linear 
system,  was  modeled  in  which  one  portion  of  the  system  repre¬ 
sented  a  unit  gain  and  the  other  portion  of  the  system 
represented  a  unit  gain  and  a  time  delay.  Three  cases  are 
analyzed  here  where  the  correlation  or  the  coherence  of  the 
inputs  has  been  varied.  The  second  type  of  experiment 
involved  a  three  mput-sirgle  output  linear  system  where  the 
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third  input  was,  in  effect,  treated  as  noise  That  is  the 
computations  are  performed  as  if  the  third  input  was  not 
known  to  exist  and  the  amount  of  its  correlation  with  the 
other  inputs  was  varied  and  its  effect  on  the  output  was 
investigated.  Three  additional  instances  are  described  here 
for  Case  II  in  which  the  gain  factor  for  the  third  system 
was  varied.  Hence,  a  total  of  six  cases  three  of  type  I 
and  three  of  type  II  were  run  for  the  computational 
experiments . 

The  next  extension  to  these  computational  experiments 
will  be  to  include  computations  of  the  output  spectra  and 
cross-spectra  directly  from  the  generated  data  rather  than 
as  a  direct  analytical  function  of  the  input  data.  This 
type  of  result  will  be  one  step  more  realistic  m  that  they 
will  contain  statistical  variability.  After  those  basic 
experiments  on  small  dimensional  cases  the  second  extension 
will  be  to  higher  dimensional  cases  to  further  evaluate  the 
more  varied  and  more  complicated  cases  ,  The  computational 
results  of  six  illustrative  cases  will  now  be  described. 

Case  I,  the  two  input-single  output  model  will  be  covered 
first,  followed  by  Case  IT  the  three  input-single  output 
model . 
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CASE  I  RESULTS 


Some  general  comments  will  first  be  made.  In  digital 
spectral  analysis,  one  can  occasionally  actually  obtain  negative 

i 

spectral  estimates  when  certain  underlying  assumptions  are  not 

perfectly  fulfilled.  This  will  never  happen  m  the  analog  case 

since  an  analog  instrument  always  filters  and  then  squares 

which  assures  that  a  positive  quantity  is  obtained.  However, 

the  smoothinq  filters  which  are  involved  in  the  digital 

analysis  occasionally  have  negative  side  lobes  which  can 

contribute  negative  power  to  the  digital  estimates.  This 

can  sometimes  result  m  coherencies  being  larger  than  1  or 

less  than  0,  even  though  m  theory  they  are  restricted  to 

the  range  0  to  1 .  Problems  such  as  these  will  often  arise 

either  where  the  spectrum  is  rapidly  changing,  ci  where  it 

is  close  to  zero.  This,-  in  fact,  has  occurred  m  some  of 

the  computational  experiments  performed  for  this  report. 

Plots  of  the  power  spectra  of  the  basic  input  processes 

G  (f)  and  G  (f)  are  illustrated  m  Figure  II1-1  The 
11  zz 

near  zero  power  m  the  higher  frequency  range  of  G  (f)  did 
actually  result  in  some  negative  spectral  estimates  and  some 
coherencies  larger  than  1  at  the  higher  frequency  points. 

The  spectra  are  plotted  at  frequency  points  of  0.2  to  1.7  cps 
and  the  bandwidth  of  the  estimates  is  .1  cps.  Hence,  these 
are  statistically  independent  spectral  points  which  are 
plotted . 
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2.1 


CASE  1(a)  RESULTS 


For  Case  I  (a),  the  parameter  a  was  taken  as  equa  •_  to 
zero  which  results  in  the  two  inputs  being  uncorrelated  or 
incoherent.  For  this  situation,  theory  indicates  that  both 
the  "true'  and  "measured"  frequency  response  function  estimates 
are  equal.  That  is,  one  can  use  either  the  correct  two- 
dimensional  theory  or  the  incorrect  one-dimensional  theory 
to  obtain  unbiased  frequency  response  function  estimates. 

This  is  illustrated  numerically  in  Figure  III-2  which  gives 
plots  cf  the  gain  and  phase  factors  for  the  true  and  measured 
versions  of  the  frequency  response  functions.  The  linear 
phase  shift  which  appears  in  the  phase  angle  plot  indicates 
the  time  delay  that  was  induced  for  the  second  input.  The 
facet  of  the  theory  that  is  not  illustrated  by  these  plots 
is  the  fact  that  the  statistical  variability  is  increased 
when  there  is  a  second  uncorrelated  input  which  is  being 
neglected.  The  statistical  accuracy  of  the  results  can  be 
shown  to  be  proportional  to  the  coherency  between  the  input 
and  the  output.  Hence,  if  one  thinks  of  the  second  input  as 
noise,  a  large  noise-to-signal  power  ratio  will  decrease  the 
coherency  between  the  first  input  and  the  output,  and  hence, 
decrease  the  statistical  accuracy  of  the  results  . 

The  coherencies  between  the  first  input  and  the  output , 
and  the  second  input  and  the  output  are  illustrated  in  Figure 
111-3,  and  reflect  the  equation  given  below  which  applies  for 
this  case  of  unit  gain  factors.  As  can  be  seen  from  this 
equation,  the  noise-to-signal  power  ratio  appears  in  the 


isE 


denominator,  and  hence  as  this  becomes  large,  the  coherency 
becomes  small. 


(f) 


1  + 


(f) 

(f) 


(1) 


The  output  v ( t )  in  this  case  is  considered  to  be  made  up  of 

a  signal  x^(t)  and  a  noise  process  x^Ct).  By  comparing 

Figure  II1-3  with  Figure  III-l,  one  sees  that  as  G  (f),  which 

z  z 

is  in  this  case  equal  to  G22(f),  becoiT>es  small,  then  the 
coherency  becomes  large.  Conversely,  when  the  coherency 
between  the  second  input  and  y(t)  is  computed,  the  reverse 
happens,  as  is  illustrated  in  the  second  part  of  Figure  III-3. 

The  ordinary  coherencies  between  the  inputs  and  the 
output  indicated  in  Figure  III-3  can,  of  course,  be  misleading 
since  there  are,  in  fact,  linear  relations  between  each 
input  and  the  output,  and  hence,  the  coherency  is  really  1. 

This  true  linear  relation  is  illustrated  by  the  plots  of  the 
partial  coherencies  given  in  Figure  III-4  which  are,  in  fact  ,  1. 
The  reason  these  go  to  1  is  that  they  are  coherencies  between 
the  input  and  the  output  after  the  best  linear  estimate  of  the 
remaining  variable  is  subtracted  out.  In  the  case  of  the 
partial  coherency  between  x^(t)  and  y(t),  the  effect  of  x^(t) 
is  subtracted  out  of  y(t)  which  means  that  the  residual  portion 
of  y(t)  will  be  directly  linearly  related  to  the  input  x^(t), 
and  this  is  illustrated  by  the  partial  coherency  between  x^(t) 
and  y(t).  The  exact  same  reasoning  applies  to  the  partial 
coherencies  between  x7(t)  and  y(t). 
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The  interpretation  of  the  partial  coherency  between 
the  two  inputs  x^(t)  and  x^(t)  becoming  unity  is  slightly 
more  involved.  The  fact  that  this  partial  coherency  becomes 
1  even  though  there  is  no  apparent  direct  linear  relation 
between  x^(t)  and  x^(t)  reflects  the  fact  that  there  is  an 
indirect  path  via  y(t)  through  which  x^(t)  and  x2(t)  are 
related.  Hence,  although  one  could  not  detect  a  relation 
between  x^  and  x^  if  one  was  measuring  only  these  two 
variables,  a  relation  between  them  does  exist  when  the  third 
variable  y  is  measured  which  is  related  to  both  of  them. 

2.2  CASE  1(b)  RESULT  3 

The  parameters  for  this  case  are  the  same  as  for  Case 
1(a)  except  tnat  now  the  inputs  are  moderately  correlated.  The 
parameter  a  which  relates  x.,  and  x^  was  chosen  to  be  0.4  which 
results  in  the  second  input  x^(t)  being  made  up  of  0.4  x, (t) 
plus  0.6  z(t).  The  major  differences  that  will  occur  here 
are  errors  in  the  computed  values  of  the  frequency  response 
function  if  just  two  signals  are  used,  that  is,  the  single 
input  and  a  single  output  and  the  effects  of  a  third  are 
ignored.  Now,  not  only  the  statistical  variability  of  the 
results  is  increased  but  also  there  are  outright  bias  errors 
that  occur  in  the  computations. 

As  can  be  seen  from  Figure  III-5,  the  differences  between 
the  true  and  measured  gain  factors  for  the  x^(t)  -  y ( t )  relation 
definitely  exist  but  I H  ( f ) |  is  not  grossly  in  error.  This 
reflects  the  fact  of  the  ignificant  but  moderate  correlation 
between  the  two  inputs.  Considerable  error  occurs  in  the 
frequency  response  function  between  the  second  input  x(t)  and 
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Frequency  Response  Functions  -.o  Case  I  (b) 


the  output  y(t)  since  in  this  case  the  output  y(t)  is  made  up 

primarily  of  x^Ct).  Hence,  large  errors  are  introduced.  This 

is  illustrated  in  both  the  gain  and  the  phase  factors  of  the 

frequency  response  between  x^(t)  and  x  (t)  as  illustrated  in 

Figure  The  ordinary  coherency  between  x.,(t)  and  x  (t)  now 

reflects  the  correlation  that  was  induced  between  the  two 

inputs.  For  Case  1(a),  it  was,  of  course,  zero  for  all 

frequencies;  however,  now  for  Case  1(b)  it  becomes  quite  large  - 

in  fact,  nearly  one  for  some  of  the  higher  frequency  ranqes . 

Actually,  values  very  slightly  larger  than  one  were  obtained 

for  the  two  highest  frequency  points  which  reflected  the 

negative  spectral  estimates  that  were  obtained  in  G  (f).  The 

zz 

ordinary  coherencies  between  each  of  the  inputs  x^(t)  and 
x^(t)  with  the  output  y(t)  exhibits  some  more  behavior  to  the 
ordinary  coherency  between  the  two  inputs.  All  of  this  is 
illustrated  in  Figure  III-6.  Figure  III-7  again  gives  the  plots  of 
the  ordinary  coherencies,  which  again  all  turn  out  to  be 
unity.  This  reflects  the  existence  of  the  true  linear  relations 
between  the  variables.  Again,  the  full  linear  relationship 
between  the  two  inputs  x^(t)  and  x^(t)  i-s  not  totally  a  direct 
relation,  but  is  partially  the  relation  induced  via  the  third 
variable  y(t) 

2.3  CASE  1(c)  RESULTS 

Case  1(c)  carries  Case  1(b)  one  step  further  in  that 
it  differs  only  in  that  the  inputs  are  more  highly  correlated. 

For  this  case  the  parameter  a  is  chosen  to  be  equal  to  0.3 
which  induces  a  stronger  correlation  between  the  two  inputs 
x^(t)  and  x^(t).  This  will  be  reflected  in  larger  errors  in 


III  -  7 


the  computation  of  the  frequency  response  function  between 
x^(t)  and  y(t).  Figure  III-8  gives  the  gain  and  phase  factors 
of  the  true  and  measured  versions  of  the  frequency  response 
functions  of  Case  1(c). 

As  can  be  seen  by  comparing  the  gain  factor  for  the  true 
and  measured  versions  of  the  frequency  response  function,  there 
is  a  considerable  difference.  As  a  matter  of  fact,  almost  an 
80%  error  is  introduced  at  the  lower  frequencies .  The  phase 
factor*  is  also  in  error  by  a  fairly  considerable  amount  going 
to  about  .88  radians  at  the  higher  frequency  end  which  would 
indicate  a  phase  change  of  almost  50°  when  done,  in  fact,  exists. 
The  frequency  response  function  between  x2(t)  and  y(t)  is  not 
in  error  in  this  case  quite  as  much  as  it  was  in  Case  1(b). 

This  is  due  to  the  fact  that  in  this  case  more  of  the  output 
y(t)  results  from  the  second  input  ^(t).  The  coherence  functions 
are  all  very  close  to  unity  as  can  be  seen  from  inspecting 
Figure  III-9.  The  coherence  between  the  two  inputs  is, of  course, 
high  because  a  very  high  correlation  was  induced  between  the 
two  inputs  and  their  construction.  The  correlation  between 
each  of  the  inputs  and  the  Oucput  is  also  high  because  in  this 
case  nearly  all  of  each  input  is  made  up  of  x, (t)  and  the  output 
is  hence  also  made  up  of  x^ft).  Hence,  y(t)  is  highly  coherent 
with  both  x^(t)  and  x2(t). 

Again,  as  in  the  previous  case,  the  partial  coherencies 
are  all  unity  since  all  the  inputs  are  accounted  for  ard 
subtracted  out  when  these  are  computed. 


Ill 
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Figure  III-10 


3.  CASE  II  RESULTS 


Case  II  differs  from  Case  I  in  that  a  third  input  is 
generated  in  the  data;  however,  the  computations  are  performed 
assuming  only  a  two  input-single  output  system.  This  will 
correspond  to  the  physical  situation  where  one  is  unknowingly 
neglecting  to  measure  a  contributing  variable.  For  this  case, 
the  correlation  between  the  three  inputs  was  maintained  constant 
in  that  the  parameters  a  and  3  were  set  to  0,4  and  0.6,  re¬ 
spectively,  for  all  three  versions  of  Case  II.  The  contribution 
of  the  extraneous  input  x^(t)  to  the  output  was  controlled  by 
varying  the  gain  factor  a^  between  three  values  which  were 
0.2,  0.5,  and  1.0.  For  small  a^,  x^(t)  can  be  thought  of  as 
correlated  measurement  noise.  For  the  larger  values  cf  a^, 
x^(t)  can  be  thought  of  more  as  a  third  input  * hich  was  not 
being  properly  considered. 

One  major  difference  that  will  occur  throughout  the 
Case  II  results  is  that  the  partial  coherencies  will  not  come 
out  to  be  perfectly  unity  since  in  this  case  not  all  the  inputs 
are  being  accounted  for.  However,  in  general,  they  will  come 
out  higher  than  the  corresponding  ordinary  coherencies  since 
they  will  more  closely  reflect  the  true  situation  in  that  as 
many  as  possible  of  the  variables  are  being  accounted  for 
properly. 

3.1  CASE  11(a)  RESULTS 

This  case  represents  small  contributions  from  the  third 
input  to  the  output.  That  is,  the  coefficient  a^  which  controls 
the  contribution  of  the  third  input  to  the  output  y(t)  is  set 
equal  to  0.2  for  this  case.  Figure  ill-ll  presents  the  true  and 
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quency  Response  Funct 


measured  frequency  response  functions  for  this  case.  Here  it 


should  be  noted  that  the  so-called  true  gain  factor  will  be  in 
error  although  not  so  much  in  general  as  the  measured  gain  factor. 

This  is  due  to  the  existence  of  the  third  input  which  is  not 
being  accounted  for.  The  true  gain  factor  for  the  frequency 
response  between  x^(t)  and  y(t)  in  this  case  actually  was  1.12, 
whicn  is  an  error  of  15%  from  the  expected  value  of  1.9.  As  can 
be  seen  from  Figure  III-ll,  the  measured  value  varies  with 
frequency  althought  not  too  greatly.  The  true  gain  factor  for 
the  second  frequency  response  function  is  computed  correctly 
as  1.0  since  the  second  input  x^(t)  is  independent  of  the 
extraneous  input  x3(t).  This  part  of  Case  II  (a)  is  therefore 
similar  to  Case  I  (a) .  The  measured  frequency  response  function 
gain  factor  is,  however,  in  error  by  a  considerable  amount  which 
is  to  be  expected  since  a,  which  dictates  the  correlation 
between  x„(<_)  and  x^(t),  is  fairly  small  and  hence  the  output 
y (t )  is  made  up  primarily  of  x  (t)  rather  than  x  (t). 

The  ordinary  coherence  functions  are  illustrated  in 
Figure  IIT-12  for  Case  II  (a) .  As  can  be  seen  they  are  consider¬ 
ably  in  error  at  the  lower  frequency  ranges,  but  approach  the  correct 
values  of  near  unity  at  the  higher  frequency  ranges.  This  is 
due  mainly  to  the  fact  that  the  power  spectrum  for  x^(t)  dominates 
at  the  higher  frequency,  and  since  the  output  is  made  up  prima¬ 
rily  of  x^(t),  the  high  coherence  come  through. 

The  partial  coherence  functions  are  shown  in  Figure  III-13 
and  as  can  be  seen  they  are  very  close  to  one  throughout  the 
entire  frequency  range.  But  they  are  definitely  slightly  less 
than  one.  This  reflects  the  fact  that  the  third  extraneous 
input  that  is  not  being  considered  has  relatively  little  effect 
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on  the  output  in  this  case.  Larger  discrepancies  will  occur 
for  Cases  11(b)  and  11(c)  which  are  discussed  subsequently. 

3.2  CASS  II  (b)  RESULTS 

Case  11(b)  differs  from  Case  11(a)  in  that  the  effect 
of  i.ne  extraneous  third  input  x^(t)  is  allowed  to  have  larger 
effect  on  the  output  y(t).  This  is  accomplished  by  setting  the 
parameter  a^  equal  to  0.5.  This  will  cause  slightly  larger 
biases  in  just  about  all  of  the  computations.  The  true  gain 
factor  will  be  in  error  by  somewhat  more  for  the  frequency 
response  function  between  x^(t)  and  y(t)  as  will  the  various 
measured  values  and  the  coherencies.  The  frequency  response 
functions  are  illustrated  in  Figure  III-14. 

In  Figure  III-14,  the  value  for  the  gam  factor  for  the 
true  frequency  response  between  x^(t)  and  y(t)  actually  comes 
out  to  1.3  instead  of  the  expected  value  1.0.  This  indicates 
the  increasing  bias  error  that  occurs  when  the  effect  of  the 
third  unaccounted  for  input  is  magnified.  However,  since  this 
third  input  is  uncorrelated  with  the  second  input  x^t),  the 
gain  factor  or  the  true  frequency  response  function  between 
x^(t)  and  y(t)  remains  correctly  computed.  This  figure  also 
illustrates  the  slightly  larger  errors  that  occur  in  the 
measured  phase  factors  and  gain  factors . 

The  coherence  functions  are  again  in  error,  of  course, 
and  are  illustrated  in  Figure  III-15  and  III-16.  The  errors  in 
the  ordinary  coherence  functions  remain  about  the  same  as  in  Case 
1(a);  however,  this  time  the  errors  in  the  partial  coherence 
functions  b  gin  to  be  increased  due  to  the  more  pronounced  effect 
of  the  extraneous  input.  In  fact,  at  the  three  highest  frequency 
points,  meaningless  results  begin  to  be  obtaired  due  to  the 
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fact  that  the  power  spectrum  of  z(t)  which  serves  to  make 
up  x^ftj  and  part  of  x^(t)  has  two  negative  estimates  which 
tend  to  invalidate  the  computations  at  that  point.  In  general, 
the  partial  coherence  functions  are  now  deviating  farther  from 
unity  than  previously  which  is  to  be  expected.  Of  course, the 
fact  that  they  are  deviating  farther  and  farther  from  unity 
can  be  used  as  an  indicator  of  the  magnitude  of  the  third 
neglected  input  which  should  be  one  of  the  general  applications 
for  the  partial  coherence  function, 

3.3  CASE  life)  RESULTS 

Case  life)  differs  from  Case  Ufa)  and  11(b)  in  that 
the  effect  of  the  third  extraneous  input  is  magnified  further 
yet.  In  this  case,  the  parameter  a^  which  controls  the  con¬ 
tributing  as  much  as  each  of  the  other  two.  Hence,  one  expects 
the  true  gain  factor  for  the  frequency  response  function  between 
x^ft)  and  y(t)  to  be  in  error  by  a  larger  amount  yet.  This, 
in  fact,  happens  as  do  the  increased  errors  in  the  partial 
coherence  functions. 

The  gam  and  phase  factors  for  the  true  and  measured 
frequency  response  functions  are  plotted  in  Figure  III-17 .  In 
this  case,  the  gain  factor  for  H^ff)  computed  to  be  1.6, 
which  is  a  60%  error  from  the  expected  value  of  1.0.  The 
measured  gain  factor  actually  varies  from  almost  2.0  to  about 
1.4,  which  is  a  considerable  error.  The  true  gam  factor 
computations  for  H  (f),  of  course,  remains  correct  at  1.0 
since  the  third  input  x,(t)  is  still  uncorrelated  with  x^Ct). 
However,  the  measured  gam  factor  varies  considerably  as  does 
the  measured  phase  factor. 
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The  ordinary  coherence  functions  and  partial  coherence 
functions  for  Case  11(c)  are  illustrated  in  Figures  HI-18  and  111-19. 
The  over-all  errors  in  the  ordinary  coherence  functions  remain 
much  the  same  as  in  the  previous  two  cases.  However ,  the  errors 
in  the  partial  coherence  functions  become  more  and  more  pro¬ 
nounced,  again  exhibiting  the  increased  effect  of  the  third 
input  x^(t).  In  fact,  in  some  cases,  the  partial  coherence 
functions  are  in  error  more  than  the  ordinary  coherence  functions 
which  is  due  to  the  relatively  large  effect  of  the  third  input. 

The  same  instability  in  the  last  two  or  three  points  of  the 
partial  coherence  function  which  results  from  negative  spectral 
density  estimates  again  occurs  for  this  case. 

This  concludes  the  presentation  of  the  numerical  results 
for  the  computational  experiments.  The  over-all  results  come  out 
essentially  exactly  as  expected .  The  Case-I  computations  illustrate 
that  the  frequency  response  functions  are  deter,  ined  correctly 
when  all  the  inputs  are  accounted  for.  In  addition,  the  Case  I 
results  illustrate  numerically  the  inaccuracies  that  occur  when 
only  a  single  input  and  a  single  output  are  used  in  computing 
frequency  response  between  an  input  and  output  function.  Case  II 
illustrates  the  degeneracy  that  occurs  when  a  third  input  is 
introduced  which  is  unaccc  ed  foe  in  the  computations.  In 
general,  the  increasing  inaccuracies  which  occur  with  an  increasing 
effect  due  to  a  third  unaccounted  for  input  are  demonstrated 
numerically . 
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CONFIDENCE  BANDS  FOR 

THREE  DIMENSIONAL  LINEAR  SYSTEM  STATISTICS 

1 .  INTRODUCTION 

The  general  form  of  the  distributions  of  frequency 
response  function  and  coherence  function  estimates  for  the 
multivariate  Gaussian  case  has  been  developed  in  Ref.'  lj  and 
in  more  detail  in  other  as  yet  unpublished  MAC  memos. 

Comprehensive  tables  for  the  exact  sampling  distribution 

of  coherence  functions  based  on  relatively  small  number  n  of 

degrees-of-freedom  (d.f.)  (up  to  n  =  25)  and  dimension  (total 

number  of  records  of  data)  up  to  p  =  10  have  been  computer 

generated,  Ref.  ^2~j.  a  study  has  been  performed  at  MAC  under 

a  USAF  concract  to  develop  Gaussian  approximation  for  sample 

coherence.  The  results  of  this  study  are  not  yet  published 

but  are  summarized  in  this  report.  A  hey  result  of  Ret.^1 

r  "I 

is  that  all  types  of  sample  coherencies  (marginal  ^ordmaryj, 
conditional  ^partialj,  and  multiple)  have  identical  distributions 
with  certain  adjustments  in  the  d.f.  which  is  employed.  Hence, 
one  set  of  tables  and  one  Gaussian  approximation  is  sufficient 
for  all  types  cf  coherence. 

The  F  distribution  arises  in  the  computation  of  the 
frequency  response  function  estimates.  This  distribution  is, 
of  course,  well  tabulated  hence  confidence  limits  are  obtained 
in  a  straightforward  manner. 
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2  .  CONFIDENCE  LIMITS  FOR  COHERENCE  FUNCTIONS 

r  1 

The  tables  of  Ref.  '  2 j  give  the  (cumulative)  distri¬ 
bution  function  (c.d.f)  as  a  function  of  n  and  p  when  n  <  25 
and  p  <  .10.  For  ordinary  coherence  between  two  variables  era 
enters  the  tables  with  p  =  2  and  n  =  2BT  as  the  number  of  d.f. 
in  the  spectrum  and  coherence  estimates.  The  quantity  B=l/r 

max 

is  the  analysis  bandwidth  and  T  is  the  record  length. 

For  the  test  cases  which  have  been  run  for  the  seismic 
experiments  the  d.  f.  are 

n  =  2(.lcps)  (250  sec.)  =  50 

The  confidence  limits  for  sample  coherence  may  not  be  obtained 
from  these  tables  since  the  limits  are  exceeded. 

A  transformation  which  is  a  very  accurate  normal 

2 

approximation  for  moderate  values  of  true  coherence  (  .4  <  y 
<  .95)  and  useful  for  more  extreme  values  is 

z  =  tanh  1  y  =  -^  tn  ^  ^ -  (1) 

where  y  is  the  positive  square  root  of  sample  coherence.  It 
has  been  empirically  determined  that  z  is  approximately  normal 
with  an  approximate  mean  value 


M 

z 


1  +  V 
1  -  y 


+ 


-P— ■ — l _ 

2  ( n  -  p  +  1 ) 


(2) 


where  y  is  the  square  root  of  true  coherence .  The  variance  of 
z  is  approximately 


2  _  _ 1_ _ 

°z  2  (n  -  p  +  1) 


(3) 
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From  the  above  equations,  if  one  is  given  a  value  of 


2 

y 


then 


p!  h  -  o  Z.  .  <  Z  <  Yi  +  o  z  , 

L  2  Z  1-0/ 2  —  z  Z  l-a/2. 


=  1 


(4) 


Let  ^  be  broken  into  two  parts 


C  =  ~  In  ]  +  Y 
2  1  -  y 


(5) 


b 


_ P  -  1 _ 

2(  n  -  p  +  1  ) 


(6) 


Then,  Eq.  (4)  can  be  rearranged  to  give  (1  -  a)  confidence 
bounds  on  y  which  are  defined  by  the  inequality 

tan  h  (z  -  b  -  Zl  _  a/2az)  <  Y  <  tan  h  (z  -  b 

(7) 

where  Z,  is  the  (  1  -  a/2)  percentile  of  the  normal  distri- 

l-a/2 

but  ion. 

The  adjustment  that  has  to  be  made  for  partial  coherence 
is  straightforward.  In  general,  one  must  subtract  from  the  d.f. 
of  the  analysis  the  number  of  variables  whose  effect  has  been 
subtracted  out.  In  general,  one  uses  n'  d.f.  where 

n'=n-(p-2)  ;  n-2BT  (8) 

In  the  two  input  single  output  case,  one  uses  n'  =  n  -  (3-2'  = 
(n  -  1  )  which  will  be  a  minor  adjustment  for  large  d.f.  Table 


I  gives  a  few  values  of  the  99%  confidence  limits  for  2BT  =  50. 
The  values  given  in  Table  I  are  plotted  in  Figure  1.  These 
apply  for  both  the  ordinary  and  partial  coherence  functions 
computed  in  the  seismic  experiments. 
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3 .  CONFIDENCE  LIMITS  FOR  FREQUENCE  RESPONSE  FUNCTIONS 


The  confidence  bands  for  the  frequency  response  functions 

depend  on  the  sample  multiple  coherence  function,  the  sample 

multiple  coherence  function  between  the  inputs  (ordinary 

coherence  for  two  inputs),  the  sample  frequency  response  functions, 

and  the  sample  input  and  output  spectral  densities.  Alternatively, 

the  confidence  bands  may  be  given  in  a  form  which  is  a  function 

of  the  sample  frequency  response  functions,  the  sample  conditioned 

spectral  density  of  the  output  y(t)  and  the  elements  of  the 

inverse  of  the  input  sample  spectral  density  matrix.  It  is  this 

second  form  which  is  most  convenient  here  and  will  be  described. 

r 

Consider  the  two  input  ,  x  (t)  and  x  (t)  single  output 
r  -)  1  2  J 

Ly(t ) j  system  being  studied.  Let  the  matrix  of  spectral 

densities  be  represented  by 


A  A  A  A 

where  S  is  2x2, S  is  2x1, S  is  1x2,  and  S  is 

xx  xy  yx  yy 

lxl.  The  function  argument  (f)  is  dropped  for  notational 

A 

simplicity.  Also,  S  is  simplified  notation  for  the  cross- 

a  12 

spectral  density  S  between  the  two  inputs  x  (t)  and  x  (t), 

Ajk  xiX2  a  _i  ^  2 

etc.  Let  SJ  represent  elements  of  S  .  The  formulas  for 

XX  XX 
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confidence  bands  may  then  be  written  in  the  form 


±po 


The  definitions  of  the  above  quantities  are  as  follows, 
the  "conditioned"  output  spectral  density  is 


AAA  A-1A 

S  |  =  s  -  s  s  s 

yjx  yy  yx  xx  xy 


A  A  A  A  A 

=  S  -  S  H  -  S  H 

yy  yi  i  y2  2 


The  term  a^is  obtained  as  follows.  First  letting  a  prime 
denote  transpose,  define  the  quantity 


A  A  A  A 

B  =  (H  -  H)  S  (H  -  H) 

XX 


where 


H  5  (H1,  H  ) 

Then  it  can  be  shown  that  the  quantity 


n~P 

F  =  - 

P 


A 

B 


A 


is  distributed  as  an  F  statistic  with  2p  and  2(n-p)  d.f. 
aQ  is  that  constant  which  satisfies 


(10) 

First 

(ID 

(12) 

(13) 

(14) 

Now 
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Hence,  one  chooses  a  confidence  level  p^.  Then  from  tables  of 

the  F  distribution  one  finds  the  percentile  value  F  2p,  2(n-p)~j 

r  i  r  n  P0  '  ' 

and  this  quantity  is  ^(n-p)/pj  j^(l/a  )-lj.  Operations  on  Eq.  (15) 

lead  one  to  Eq  (10)  which  provides  the  simultaneous  confidence 

levels  given  as  radii  around  the  points  H  and  in  the  complex 

plane . 


The  confidence  formula  given  by  Eq.  (10)  may  be 
interpreted  in  terms  of  the  gain  factor  H(f)  and  the  phase 
factor  $(f)  by  considering  the  sketch  in  Figure  2. 


Confidence  Bands  for  Gain  and  Phase  Factors 

Figure  2. 


From  Figure  2  one  sees  that  |H  |  is  bounded  by 

,  a  ,  A  A  |  a  .  i 

jH^|+  r^  where  r^  =  |H  -  H, |  and  the  phase  factor  is  bounded 

with  confidence  by  where 


A$  =  arc  sin  — — 

1  Ih. 


(16) 


All  these  statements  hold  simc.1  aneously  for  Ih^I  and  |#2). 

The  preceding  confidence  formulas  may  be  directly  extended 
to  p  variables  by  defining  the  matrix  S  of  Eq.  (9)  for  p  vari¬ 


ables  . 


a  a  a 

S  S  •  •  •  s 

11  12  l,p-l 


s  =  * 


P-1,1 


P-1, P-1 


Y'P-1 


(17) 


p-i.  y 


Equation  (10)  becomes 


-  K'*2 


c  c11 

S  |  s 
y  I  x  xx 


I2  <(-  - 

’-1  —  V  an 


i)  i  ,  s?-1'  P-1 

/  y I x  xx 


(18) 
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and  Eq.  (11) 


written 


AAA  A -1  A 
S  |  =  S  -  S  S  S 

y|x  yy  yx  xj  xy 

AAA  A  A 

-  S  -  S  H  -  ...  -  S  H  .  . 

yy  yl  1  y,p-l  p-1  >19) 

Due  to  the  fact  that  several  variables  exist  in  the  confidence 
formulas,  it  is  not  practical  to  prepare  tables.  The  computations 
should  be  added  to  the  basic  computing  program  and  given  as  a 
part  of  the  normal  output.  A  table  of  F  values  for  varying  n 
and  p  could  be  included  in  the  computer. 
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